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Abstract 



One of the most important subjects of debate in the formation of the solar 
system is the origin of Earth's water. Comets have long been considered as the 
most likely source of the delivery of water to Earth. However, elemental and 
isotopic arguments suggest a very small contribution from these objects. Other 
sources have also been proposed, among which, local adsorption of water vapor 
onto dust grains in the primordial nebula and delivery through planetesimals 
and planetary embryos have become more prominent. However, no sole source of 
water provides a satisfactory explanation for Earth's water as a whole. In view of 
that, using numerical simulations, we have developed a compound model incorpo- 
rating both the principal endogenous and exogenous theories, and investigating 
their implications for terrestrial planet formation and water-delivery. Comets are 
also considered in the final analysis, as it is likely that at least some of Earth's 
water has cometary origin. We analyze our results comparing two different water 
distribution models, and complement our study using D/H ratio, finding possible 
relative contributions from each source, focusing on planets formed in the habit- 
able zone. We find that the compound model play an important role by showing 
more advantage in the amount and time of water-delivery in Earth-like planets. 

Subject headings: Solar system: formation, Earth, planets and satellites: formation, 
astrobiology, methods: N-body simulations 
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1. Introduction 

The origin of Earth's water is one of the most outstanding questions in planetary science. 
The locations of the regions within the early solar system where water-carrying objects 
might have contributed to the accretion of Earth is a matter of intense debate. It is widely 
accepted that the protosolar nebula was hotter and denser towards its center, and cooler 
and less dense farther out (Encrenaz, 2006). This temperature gradient significantly affected 
nebula's chemical composition and the distribution of its water and icy materials. Close 
to the center only metal and silicates condensed whereas more volatile materials condensed 
farther out from the Sun. The first solid particles had microscopic sizes. These objects 
gradually stuck to other particles and grew to larger dust grains. In an oxygen-rich region, 
debris formed carbonaceous chondrites that carry up to 10 percent water (Morbidelli et al. 
2000). Beyond the giant planets, water condensed in large quantities and formed comets, 
which are up to 80 percent made of ice (Jessberger, Kissel & Rahe, 1989). . 

Compared with these icy materials, Earth contains little amount of water. Only about 
0.02 percent of Earth's mass is in its oceans, and somewhat more water sits beneath its 
surface. An estimate of the amount of water inside Earth points to values ranging from 10® 
(O® = mass of Earth's oceans = 1.4 x 10 24 g) to 50 O® with ~ 10 O® being the value that 
is generally considered as the amount of water inside Earth's mantle (Drake & Campins, 
2006). A recent estimate of the amount of water in Earth's interior by Marty (2012) also 
agrees with these values. As suggested by this author, Earth contains ~ 4 — 120© in its 
interior. The important fact is that, irrespective of these values, Earth still has substantially 
more water than expected from a body at 1 AU from the Sun. 

In the past decade many attempts were made to explain the origin of Earth's water. 
Suggested sources can be divided into endogenous and exogenous, including water adsorbed 
by dry silicate grains in the nebula (Stimpfl, Lauretta & Drake 2004), delivery through as- 
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teroids, comets, planetary embryos, and planetesimals (Morbidelli et al. 2000; Raymond, 
Quinn & Lunine 2004; O'Brien, Morbidelli & Levison 2006; Raymond et al. 2004, 2006, 
2009; Lunine 2003; Drake & Campins 2006), and water production through oxidation of an 
hydrogen-rich atmosphere (Ikoma & Genda, 2006). These different sources may be distin- 
guished by their isotopic differences such as their deuterium-to-hydrogen (D/H) ratios as 
well as their predictions of the delivered amount of water relative to the total mass. 

Among these possibilities, comets have long been considered as an attractive exogenous 
source of Earth's water. These objects that are made of ice, are believed to likely retain 
the isotopic composition that they acquired during their formation (Drouart et al. 1999). 
Whether comets were the sole source of the delivery of water to Earth is, however, uncertain. 
The measurements of the D/H ratio of water in eight Oort Cloud comets are on average twice 
that of the Standard Mean Ocean Water (SMOW) and fifteen times the value of the D/H 
ratio in the early solar nebula (Table 1). Although the original value of the D/H ratio of 
Earth's water is unknown, and it is unclear how that value changed during the geophysical 
and geochemical evolution of Earth (for instance, as shown by Genda & Ikoma (2008), the 
Earth's D/H ratio could have increased by a factor of 9 had Earth had a massive hydrogen 
atmosphere that experienced slow hydrodynamic scape), many researchers have used the 
incompatibility between the D/H ratio in comets and in SMOW as an argument to rule out 
comets as the main source of the delivery of water to Earth. 

To resolve this issue, Owen & Bar-Nun (1995) and Delsemme (1998) suggested that 
comets formed in Jupiter's region may be the source of Earth's water. These comets have 
lower D/H ratios since, compared to the comets from the Oort Cloud, they have formed 
in a warmer environment [e.g., the D/H ratio of the Jupiter-family comet 103P/Hartley 2 
is almost the same as that of the SMOW, (Hartogh et al. 2011)]. However, models of the 
dynamical evolution of solar system do not seem to support this idea. As shown by Morbidelli 
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et al. (2000), assuming Jupiter and Saturn were formed in their current orbits, cometary 
material delivered to Earth might have been originated from the regions beyond the orbit of 
Uranus as the comets formed in the Jupiter region would have very short lifetimes implying 
that their probability of their collision with Earth has been very low. Gomes et al. (2005) 
who showed that a cometary influx to Earth based on the dynamics of the outer planets 
delivers no more than ~ 1.8 x 10 23 g of material which is only about 6 percent of the current 
mass of Earth's oceans. As shown by Morbidelli et al. (2000) and Dauphas et al. (2000), 
comets constitute a minor source of Earth's water and the contribution of cometary water 
is smaller than ~ 10%. It is, however, important to note that as shown by Ipatov & Mather 
(2007), if the effect of the dynamics of the outer planets on the scattering of Jupiter family 
comets are studied during the times that are earlier than the time considered by Gomes et 
al. (2005), the probability of the collision of these objects with Earth would be large enough 
to justify the delivery of all Earth's ocean entirely through these bodies during formation of 
the giant planets. 

Another challenge to the notion of water-delivery by comets comes from the analysis of 
noble gases and other elemental isotopic ratios. For instance, Swindle & Kring (2001) have 
argued that based on their analysis of noble gas ratios, comets could not have been able 
to supply a significant fraction of Earth's water unless either the comet delivery of water 
occurred in the first 100 Myr of Earth's history, or water was delivered by comets from 
regions other than the Oort cloud. Based on elemental and isotopic arguments, Drake & 
Righter (2002) were able to confirm the findings of Swindle & Kring (2001) and limited the 
contribution of comets to Earth's water to 50 percent. 

Dynamical simulations and other isotopic analyses suggest that a more reliable value 
for the contribution of cometary water to Earth is probably smaller than 10-15 percent 
(Morbidelli et al. 2000; Drake & Righter 2002). This is in agreement with the finding of 
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Hutsemekers et al. (2009) who used 15 N/ U N isotopic ratio and showed that the amount of 
cometary water delivered to Earth could not be more than 9%. It is, however, necessary to 
caution that studies based on isotopic analysis can lead to discussions as to which measured 
value of the D/H ratio in comets is a true representative of their bulk composition since the 
primary source of information about comets continues to be the studies of their comae. As 
shown by Weirich, Brown & Lauretta (2004), the D/H ratio would be expected to rise during 
diffusion and sublimation. This was also shown in an experiment by Schmidt et al. (2005) 
where these authors studied comet sublimation and showed that an upward trend is observed 
in the D/H ratio in the evaporated material independent of their bulk composition. These 
authors argued that in order to obtain a more realistic value of the D/H ratio in comets, it 
is important to understand the amount of bulk ice in a comet's coma. 

Another possible exogenous source of Earth's water is the primitive asteroids. The 
D/H ratios of the individual carbonaceous chondrites range from 1.2 xl0~ 4 to 3.2 xl0~ 4 
(Lecuyer et al. 1998), a value that is very close to the D/H of SMOW, implying that hydrated 
carbonaceous asteroids originated from the primordial asteroid belt might have been able 
to provide terrestrial planets' water. Water-carrying asteroids from the outer asteroid belt 
(> 2.5 AU) can be scattered into the terrestrial region due to their interactions with giant 
planets and planetary embryos. While this presents a viable mechanism for the delivery of 
water, models of the dynamical evolution and sculpting of asteroid belt indicate that its 
efficiency is low and its contribution to the total water budget of Earth may be very small. 
The original simulations of Morbidelli et al. (2000) in which planetesimals were considered 
to be massless particles, for instance, suggest only 0.13% of primitive asteroids would have 
been accreted by Earth (Morbidelli et al. 2000). Assuming 10% water content for a typical 
hydrated carbonaceous asteroid, and no loss of mass and water in each collision, this rate 
of accretion requires the mass of the asteroid belt beyond 2.5 AU to be at least 4 times 
the mass of Earth in order for primitive asteroids to deliver the current amount of water to 
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Earth (Morbidelli et al 2000). More advanced, high- resolution simulations by Raymond et 
al. (2007) show that when the planetesimals are considered to have mass, the contribution 
of primitive asteroids could be as high as ~ 5% of the initial population. 

In these preliminary simulations by Morbidelli et al (2000), a second issue with primitive 
asteroids as the source of Earth's water is the time of water-delivery. While the mass- 
requirement of water-delivery through primitive asteroids is comparable with the estimates 
of the original mass of the asteroid belt, the time of the delivery is not comparable with the 
time of terrestrial planet formation. For instance, considering primitive asteroids as massless 
particles, Morbidelli et al (2000) showed that the process of the delivery of the current amount 
of Earth's water to the Earth's accretion zone will take no longer than 40 Myr. At this time, 
Earth is young with a mass < 60% of its current value, and still forming. As Earth continues 
its growth, it is subject to numerous impacts by planetesimals and several major collisions 
with planetary embryos. It is not certain (and in fact it is even doubtful) if Earth could 
have retained water that it received through the impacts of primitive asteroids (Genda & 
Abe 2005, Canup & Pierazzo 2006). In the framework of this scenario, it is important to 
note that in more recent experiments by Raymond et al. (2006, 2007, 2009) and O'Brien 
et al. (2006), where planetesimals have mass, a significant part of the delivery of asteroidal 
material to planets around 1 AU happens when these bodies have accreted more than 50% 
of their final masses. 

The requirement that the asteroid belt must have been much more massive in the 
past in order to facilitate the delivery of water by asteroids, combined with the fact that 
asteroids/planetesimals and planetary embryos (moon- to Mars-sized bodies formed shortly 
after the formation of planetesimals) are scattered into varieties of orbits as a result of 
interacting with one another and with giant planets suggests that the delivery of water 
to Earth must have occurred during the entire course of Earth's formation through the 
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impacts of hydrated planetesimals and water-carrying planetary embryos. These objects 
originated from the outer asteroid belt and deposited their water contents when impacting 
Earth (Morbidelli et al 2000; Raymond et al, 2004, 2006, 2007, 2009; O'Brien et al, 2006). 
This model is widely accepted as the main mechanism for the delivery of water to Earth and 
points to the outer asteroid belt as the reservoir of the water-carrying materials. 

The main drawback for considering asteroidal water as a source of Earth's water is 
that the Earth's primitive upper mantle has a significantly higher 187 0s/ 188 0s ratio than 
carbonaceous chondrite. In fact, Os isotopic composition of Earth's primitive upper mantle 
matches those obtained from anhydrous ordinary chondrites, and is distinctly higher than 
anhydrous enstatite chondrites (Drake, 2005). 

It is important to note that the above-mentioned models of water-delivery are based 
on the assumption that giant planets maintain their (current) orbits during the formation 
of terrestrial planets. Recently Walsh et al (2011) have shown that the delivery of water 
through hydrated asteroid and planetary embryos can also occur during the migration of 
giant planets. Known as the Grand Tack model, these authors proposed an early inward- 
then-outward migration of Jupiter and Saturn in a gas-rich phase, and showed that unlike 
the previously mentioned models, the water-delivery to terrestrial planet zone occurs when 
both planets are migrating outward and scatter water-rich bodies initially orbiting beyond 
the orbit of Jupiter, to the inner regions. 

In addition to the exogenous sources explained above, it has been suggested that the 
Earth's water might have an endogenous origin: water could have come directly from the solar 
nebula. Stimpfl et al. (2004) examined the role of physisorption by modeling the adsorption 
of water on dust grains at 1000 K, 700 K, and 500 K using Monte Carlo simulations and 
showed that grains accreted to form Earth could have adsorbed 1-3 Earth oceans of water 
prior to their accretion. As shown by these authors, the efficiency of water adsorption 
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increases as the temperature decreases. 

Although the observations of water vapor and atomic hydrogen at 1 AU around young 
stellar objects such as MWC 480 (Eisner, 2007), and the presence of forsteritic olivine in the 
dust clouds surrounding young solar-type stars are supportive of the water adsorption model 
(Muralidharan et al. 2008; Stimpfl et al. 2006; Eisner, 2007; Bethell, & Bergin, 2009), this 
scenario suffers from issues related to the retention of water during the growth of dust grains 
to larger objects. Another issue with this model is the value of the D/H ratio since the D/H 
ratio of the nebula is 7 times smaller than that of SMOW. That means, if the adsorption 
process had provided all the water on Earth, a mechanism would be required to raise the 
D/H ratio from its solar value to that of the SMOW (Drake, 2005). 

It is necessary to note that the compositional gradient of planetesimals and planetary 
embryos in the solar nebula, in particular the concentration of their volatiles depends on their 
thermal evolution during their growth. Several factors contribute to the thermal evolution of 
a body including local heating by the Sun, concentration of radioactive nuclides, proximity 
to the magnetic field of the Sun, and the final size of the body. These processes are very 
complex and their study is beyond the scope of this paper. We refer the reader to Nuth 
(2008) for more details. 

While, as explained above, the origin of Earth's water is still unresolved, the uncertain- 
ties regarding the amount of water in Earth's mantle have turn this issue into a topic of 
debate as well. For instance, while Hirschmann et al. (2005) argue for a highly wet mantle 
that contains approximately 20O ffi , the model by Smyth et al. (2006) suggests an almost 
dry mantle with no more than 2 Earth's ocean. A comprehensive model of the formation 
of Earth and origin of its water has to address these issues as well. As mentioned above, 
planetesimals and planetary embryos have been able to provide a large portion of Earth's wa- 
ter. However, other sources such as comets, primitive asteroids, and water adsorbed on dust 
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grains have also had their shares and contributions. What separates these sources from one 
another is the time of their operation. As indicated by geological data, each of these events 
occurred at a different time during the formation of Earth. In order to be able to address the 
problem of the origin of Earth's water properly, it is necessary to determine how much water 
was delivered by each of these mechanisms, and at what time during the terrestrial planet 
formation. This paper presents an attempt to address these questions. Considering that 
the contribution from comets is < 10% of the Earth's water budget, we focus our study on 
water-delivery through planetesimals and planetary embryos, and combine that with water 
adsorption on dust grains as a source of Earth's water. Our goal is to answer the following 
question; Assuming water adsorption on dust grain as a viable mechanism for the delivery 
of some of Earth's water, would the inclusion of this mechanism in the currently accepted 
model of water-delivery through planetesimals and planetary embryos be consistent with 
geological constraints on the amount of water in Earth? A positive or negative result will 
enable us to assess the validity of the scenario in which water adsorbed on the surface of 
dust grains provides the source of Earth's water. 

In the following section we discuss our model. Section 3 presents the details of our 
numerical experiment, and in section 4, we analyze the results. Section 5 concludes this 
study by presenting a summary and discussing the implications of its results. 

2. Model 

Since we are interested in studying the contribution of water from endogenous sources, 
we consider two disk models "A" and "B". The disk in model A is bi-modal. It consists of 
an inner part from 0.5 AU to 2.5 AU where it is populated only by planetary embryos and 
an outer region extending from 2.5 AU to 4 AU where it contains only planetesimals. The 
surface density of the disk is assumed to have the following two-tiered radial profile 



' S^r/IAU)- 3 / 2 , 0.5AU < r < 2.5AU 

E A (r) = < 

k S 2 (r/5AU)- 3 / 2 , 2.5AU < r < 4.0AU. 



(1) 



Model "B" is a disk of planetary embryos only, and extends from 0.5 AU to 4 AU. The 
surface density of this disk is given by 

S B (r) = E^r/IAU)- 3 / 2 , 0.5AU < r < 4.0AU (2) 

In both models, planetary embryos are distributed randomly with a mutual separation of 
5-10 Hill radii. At each semimajor axis, the mass assigned to an embryo is given by 

3/2 

(3) 

where a is the semimajor axis of the embryo, E = Sa(o), S B (a), is the disk's surface density 
at the position of the embryo, and M G represents the mass of the Sun. 

The planetesimals in the outer region of the disk model A are assumed to have a mass 
of 0.01 Earth-masses. We distributed 200 of these planetesimals between 2.5 and 4 AU 
following the surface density of the disk as given by equation (1). All planetesimals and 
planetary embryos were initially in circular orbits. The orbital inclinations of these objects 
were chosen randomly from the range of 0.0001 to 0.001 degrees. The number of planetary 
embryos in disk model A is between 60 and 80 and in disk model B is between 80 - 100. 

Our models also include giant planets. As shown by many authors (Chambers & Wether- 
ill 1998; Chambers 2001; Raymond et al. 2004, 2006, 2009), these objects play an important 
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role in the radial mixing of the disk material, and the final assembly and water contents 
of terrestrial planets. We consider two models, one with Jupiter, and one with Jupiter 
and Saturn, and assume that the orbital elements of these objects, at the beginning of our 
simulations, are similar to their current values. 

We consider a biphasic scenario for the distribution of water in both our disk models. 
We consider planetesimals interior to 0.7 AU to be dry and assume that in the region of 
0.7-2.5 AU, water is distributed according to the (endogenous) model proposed by Stimpfl et 
al. (2004). In this model, the amount of water varies according to the temperature gradient 
profile of the solar nebula. An examination of the relation between temperature and distance 
in the solar nebula as given by Clark (1998), and the relation between distance and water 
adsorbed by planetesimals and planetary embryos in the inner Solar System as given by 
(Stimpfl et al. 2004) suggests that the amount of water adsorbed by an object in the inner 
part of the disk is given by 

W = 10a - 7. (4) 

In this equation, W is the water content adsorbed (in ocean water) by terrestrial masses and 
a is the object's semimajor axis in AU. For the outer region of the disk (a > 2.5 AU) in 
both models, we use a water content representative of carbonaceous chondrites, equivalent 
to 5 percent of the mass of each object (Raymond et al. 2004, 2005a, 2005b, 2006, 2009). 
Following equation (4), a body interior to 0.7 AU is dry, at 1 AU it will have 0.07% water 
by mass, at 2 AU its water-mass fraction will be 0.3%, and at 2.5 AU, 0.42% of its mass will 
be water. This range of water-mass fraction indicates that the distribution of water in our 
model is consistent with the amount of water carried by S-type asteroids. These objects are 
the dominant bodies in the inner asteroid belt and contain ~0.1% water by mass (Abe et 
al., 2000). Despite this agreement, it is still hard to explain the presence of more anhydrous 
primordial bodies in the inner solar system if the adsorption of water vapor on dust grains 
were an efficient process in the solar nebula (Drake, 2005). 



3. Numerical Simulations 



We integrated our systems for different combinations of the disk surface density profile 
and giant planets configurations. In both disk models (A and B), we considered three 
different values for the surface density at I AU (£1): 6 g/cm 2 corresponding to the surface 
density of the inner solar system in the minimum mass solar nebula (Weidenschilling 1977; 
Hayashi 1981; Raymond et al. 2004), 10 g/cm 2 as the commonly considered maximum 
surface density at 1 AU, and 8 g/cm 2 as an intermediate value. For S 2 we only considered 
3 g/cm 2 . 

We performed 54 simulations; 36 with the disk model A and 18 with the disk model B. 
In the simulations with model A, 18 were carried out with only one giant planet (Jupiter) and 
18 included both Jupiter and Saturn. In these simulations, each disk with a different value of 
Si was integrated 6 times with slightly different initial configurations for planetesimals and 
planetary embryos. In the simulations with model B, 9 simulations included only Jupiter 
and 9 included both giant planets. In this case, we performed three simulations for each 
value of Si each with slightly different initial configurations. 

Simulations were carried out for 300 Myrs using the hybrid integrator in the N-body 
integration package Mercury (Chambers 1999). The time step of simulations was set to 5 
days. We assumed all collisions to be perfectly inelastic and conserve linear momentum. 

Figures 1 and 2 show the final results for all 54 simulations. The labels on the vertical 
axes represent the number of the simulation, the giant planet (s), and the values of Ei. The 
color of each object corresponds to the value of its water- mass fraction based on the color 
coding scheme at the bottom of each figure. The orbital eccentricity of each object is shown 
by a horizontal line corresponding to the range of the variations in the heliocentric distance of 
the body from its perihelion to aphelion. Table Al in the appendix shows the corresponding 
values of the mass, semimajor axis, orbital eccentricity, and water content of each body. As 
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a point of comparison, the inner planets of the solar system and their water contents are 
also shown. The gray area shows the habitable zone of the Sun from the model by Kasting 
et al. (1993). 

An inspection of the results show a great diversity in the mass, water content, and 
orbital configurations of the final objects. Figure 3 shows the initial and final values of the 
total masses of the planetesimals and planetary embryos in each disk model. As shown here, 
in all simulations, the disk lost a great portion of its mass. This is an expected result that is 
due to a combination of ejection (i.e., reaching a heliocentric distance of 100 AU or larger), 
and collision with Sun or the giant planets. Tables 2 and 3 show the fraction of the mass lost 
by each mechanism. As shown here, most of the mass of the disk was lost as a result of the 
ejection of objects out of the system. The ejected mass is larger in systems with two giant 
planets. Almost always more than 50% of the initial mass of the disk was ejected in these 
systems. This is due to the fact that in such systems, in addition to the direct gravitational 
perturbation of giant planets, secular resonances such as u Q excite the orbits of planetesimals 
and planetary embryos and cause many of them to be scattered out of the system. We refer 
the reader to Levison & Agnor (2003), Raymond et al. (2009) and Haghighipour et al. (2013) 
for an analysis of the effect of z/ 6 on the scattering of materials from a protoplanetary disk 
and the onset of terrestrial planet formation. At the end of our simulations, the mean change 
in the semimajor axes of Jupiter and Saturn were 5 aj ~ 0.03 — 0.05 and 5 as ~ 0.01 — 0.03, 
respectively. These values are consistent with those reported by Raymond et al. (2004). 

The strong perturbations of Jupiter and Saturn combined with the effect of secular 
resonances in ejecting material from the disk in simulations where both these two planets 
were included results in smaller number of collisions in these systems compared to those 
simulations that included only Jupiter. For instance, in the simulations with disk model 
A, the total massloss due to collision with Jupiter is ~ 0.2M e whereas in the simulations 
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with the same disk when both Jupiter and Saturn are included, this value drops to 0.06M®. 
Same results have been obtained for simulations with disk model B. The situation is, however, 
different for the collisions with the Sun. In this case, the amount of the mass lost due to the 
collision with the Sun is higher in systems that include Jupiter and Saturn. 

An interesting result of our simulations is the connection between the values of the mass, 
semimajor axis and orbital eccentricity of the final objects, and the disk surface density profile 
and giant planet (s) configurations. Figures 4 and 5 show the results in a mass- semimajor 
axis diagram. In these figures, each box represents the final results of all 6 simulations 
that correspond to its mentioned value of the surface density at 1 AU and giant planet's 
orbit. As shown here, larger planets are formed in the region between 0.5 AU and 1.5 AU 
and their mass distribution peaks at approximately 0.7-0.8 AU. This is an expected results 
that attributes to the fact that in the inner regions, the perturbing effect of giant planets 
is negligible and the dynamics of planetesimals and planetary embryos is primarily driven 
by the interactions of these bodies with one another. A consequence of the latter is that in 
a disk with a larger surface density (where more material is available in its inner region), 
the collision and accretion of planetesimals and planetary embryos result in a few but larger 
objects (Raymond et al. 2004, 2005b). In the outer region (> 1.5 AU), as shown by these 
figures, the perturbation of giant planets causes the disk material to be scattered out of the 
system and as a result, the final objects formed in these regions carry less mass and are 
smaller. This can, for instance, be seen in the simulations of model B with X = 6g/cm 2 
where because of the effects of Jupiter and Saturn, only a few objects smaller than Earth 
survived the integrations. 
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3.1. Water content 

Figures 4 and 5 also show the water-mass fraction of the final objects. It is important 
to note that the final water content of a body, in addition to the contributions from water- 
carrying planetesimals and planetary embryos initially orbiting past 2.5 AU (hereafter refer 
to as asteroidal water), depends also on the amount of the water that was initially adsorbed in 
the inner regions of disk as Eq. 4 (simulations of model B included only planetary embryos). 
Results of the simulations of the two sets of model A (Tables 4 and 5) indicate that the 
percentage of the asteroidal water received by planets formed around 1 AU varies between 
0% and 74%. Despite the large amount of water that was initially available in the disk, 
some of the final objects in this model received different amounts of water even as small as 
10®, the minimum value of Earth's water content. This is due to the fact that many water- 
carrying bodies, especially those in the outer region of the disk, were dynamically excited by 
giant planets and were scattered out of the system. 

An interesting result obtained from the simulations of model A is that despite the 
different amounts of asteroidal water and mass that were received by final planets, the 
percentage of the asteroidal water in these objects in simulations with one or two giant 
planets is the same. As shown by Table Al, the amount of asteroidal mass received by most 
of these planets is smaller than 10%. Figure 4 shows that similar to the mass-distribution of 
the final objects, the planets with more water are within 0.5-1.5 AU with their water-mass 
ratios peaking at 0.7-0.8 AU. 

Table Al also shows small planets with almost no asteroidal water implying that these 
objects accreted only a few or no planetesimals that were orbiting initially past 2.5 AU. 
Naturally, these objects would be expected to have formed in the inner region of the disk 
where the amount of asteroidal water is small or zero. However, our simulations show 
that objets without asteroidal water could also form at semimajor axes larger than 2.5 AU. 
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Probably these bodies were scattered to the outer region of the disk by interacting with 
embryos from the inner part remaining with no contribution of water from asteroidal bodies. 
See Haghighipour & Scott (2012) for the possibility of scattering of planetesimals from inner 
orbits to outer regions. 

Note that the cases 24A-JS and 13A-J (Table Al) have indeed some amount of water. 
However, the amount of water in these bodies is smaller than 0.05 Earth's oceans. Because 
of such a small amount of water, and because in the table, we only considered one digit 
after the decimal point in our entries, the water contents of these bodies do not show. The 
planet in the system 3B-J (Table Bl), on the other hand, is actually dry. That is because 
this object was initially placed between 0.5 AU and 0.7 AU which is a dry region in our 
compound model. 

3.2. Water in the Habitable Zone 

We followed Kasting et al. (1993) and considered a region for the habitable zone (HZ) 
of the Sun extending from 0.9 AU to 1.37 AU. In all our 54 simulations, a total of 39 planets 
formed inside the HZ. Tables 4, 5, and 6 show the final results. One can see from these 
tables that the masses of these planets vary between ~ 0.2M e and ~ 2.2M e with a water 
content ranging from 10 e to 380 e . The maximum water content of the planets formed in 
the HZ did not exceed the maximum value of 50O e expected for Earth (Abe et al. 2000). In 
the following, we present a detailed discussion on the amount of water carried by the final 
planets, the time of the delivery of that amount of water, and use the D/H to determine 
the amount of cometary material that needs to be delivered to the final planets in order for 
the D/H of their water contents to match that of Earth. To portray the effect of the initial 
water adsorption, we also make a comparison between the results of our model and those 
adopting the model used by Raymond et al. (2004; 2006; 2007; 2009). 



3.2.1. Water content of the final planets 



When including only Jupiter, most planets in the simulations of models A and B were 
able to carry at least 5 Earth oceans of water (Tables 4-6). However, as shown in Table 
7, the efficiency of the delivery of more than 10 Earth oceans, which by some estimates is 
the total amount of water carried by Earth (Drake et al. 2005), dropped significantly in 
our simulations when we considered the water distribution model of Raymond et al. (2004, 
2006, 2007; 2009). In the latter case, only 2 of the 13 planets (< 20%) that were formed in 
the HZ received at least 10 Earth's oceans. As a point of comparison, our compound model 
(Table 4) which included water adsorption as well, was able to deliver 10 Earth's oceans in 
11 of the 13 planets formed (> 80%) in HZ. 

A comparison between the results in Tables 4 and 5 shows that in the simulations of set 
A, the inclusion of Saturn significantly reduced the total amount of water in the final planets 
that formed in the habitable zone. However, this reduce in efficiency was still less than that 
in models in which the water was delivered only through collisions of planetesimals and 
protoplanetary bodies (Raymond's model), and did not include adsorbed water. As shown 
in Table 7, in our compound model, even when Saturn is included, 12 out of 14 planets in 
the HZ received more than 5 Earth's oceans. When considering the model without water 
adsorption (in our same simulations), this number dropps to only 2 out of the same 14 
planets. This reduction in efficiency in both models is due to effect of secular resonance 
with Saturn. This resonance increases the eccentricity of small water carrying bodies in the 
region around 2.2 AU, and causes many of these objects to be ejected from the system or 
collide with the Sun in a short time. 

Interesting results were obtained in simulations of model B. In these simulations, despite 
that our model included water adsorption, no difference was observed in the amount of the 
water carried by the final planets when Saturn was included in the simulations. The delivery 
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of at least 5 Earth's oceans was observed in roughly 50% of planets that formed in the HZ. 
This could be attributed to the fact that in these simulations, the disk model B contains 
only planetary embryos which, unlike planetesimals, are not easily excited by Saturn's u 6 
secular resonance [the embryo-embryo interactions around 2.0-2.5 AU in these simulations 
neutralize the effect of the u 6 resonance. See Levison & Agnor (2003) and Haghighipour et 
al. (2012) for more details]. Since the embryos past 2.5 AU carry 5% water, a single collision 
with a proto-planet in HZ can deliver substantial amount of water to a planet in that region. 



3.2.2. Time of water delivery 

The time of water-delivery also showed dependence on the giant planets configuration. 
For instance, in the simulations of the disk model A considering our compound model (model 
A- J, see Table 7), it took on average 19 and 36 Myr to deliver 2 and 5 Earth's oceans to the 
final planets in the HZ (Tables 4 and 7), respectively. This is while in our same simulations 
using a model of water distribution as in Raymond et al (2004, 2006, 2007; 2009), where 
water adsorption is not included, this time is approximately 26-117 Myr. When Saturn was 
included (model A-JS, see Table 7), as expected the time of delivery was prolonged. In our 
simulations (Tables 5 and 7), the timescale required for a planet in HZ to receive 2-5 Earth's 
oceans was on average 40-56 Myr using our compound model. In models without water 
adsorption, this time extends to 52-181 Myr. These results indicate that in general, our 
model is more efficient in delivering water during the entire course of the planets' accretion. 
For instance, planets formed inside the HZ from set A suffered their last giant collision on 
average after 150 Myr with a range varying from 30 Myr to 290 Myr (Table 7). 

While in the model without water adsorption, only 2 planets received 10O e on average 
within 90 Myr, and no planet received more than 150 e , in our compound model, 13 plan- 
ets received 10O ffi in about 100 Myr and 4 planets received 150 e on average in 146 Myr. 



-20 - 



Evidently, our compound model is more efficient in delivering water to many planets during 
the formation of these objects. We recall that the timescale for the formation of Earth is 
between 50 Myr and 150 Myr (Touboul et al. 2007). 

Similar results have been obtained in our simulations of disk model B (Tables 6 and 7). 
For instance, simulation of 12B-JS in Table 6 produced a planet with 2.26 Earth's masses 
carrying a total 38 Earth's oceans. From this amount of water, 17 Earth's oceans were 
strictly due to local adsorption of water in the primordial solar nebula. As shown in the 
table, water-delivery through planetesimals and planetary embryos stopped after 31.2 Myr. 

In regard to the time of the delivery of water to Earth, it is important to determine what 
fraction of the delivered water could have been retained by the Earth if the water-delivery 
occurred very early on when the Earth was only half of its current mass (Morbidelli et al. 
2000). To answer this question, we first compare the time of water-delivery for the delivery of 
first two oceans (< 20 ffi ) for both models of water distribution (i.e., with and without water 
adsorption). For the simulations of the disk model A where only Jupiter was included, the 
delivery of the first 2 oceans occurred for both models of water distribution between 10 and 
35 Myr in 50% of planets formed in the HZ. As shown in Table 7, this timescale is between 
0.5 and 62.1 Myr for our water distribution model and between 2.7 and 62.1 Myr for models 
in which water adsorption is not included. The close similarity between these timescales 
indicates that the time of the delivery of the first two oceans in both models is statistically 
indistinguishable. In fact, in several of our simulations, the timescale for the delivery of the 
first 2 Earth's oceans to planets in the HZ was equal in both water distribution models. This 
result implies that our compound model does not play a significant role in the early delivery 
of water. Given the statistical indistinguishability of the results of both models, we conclude 
that both water distribution models can result in similar loss of water in the early accretion 
of planets in the HZ. However, as the accretion of these objects proceeds, the compound 
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model gains an advantage in the amount and time of water-delivery compared to models 
of water distribution as that in Raymond et al (2004, 2006, 2007, 2009) and O'brien et al 
(2006). 

Our compound model is also capable of delivering much more water to many more 
planets in a timescale consistent with the mean-time of the last giant collision. Table 7 
shows the results. As shown in this table, for several planets in both models, water accretion 
reaches its maximum after these planets grew to 60% of their final masses. The mean-time 
of 60% accretion and the range of the time of the 60% accretion are shown in the table as 
well. In some cases, more than 50% of the final water content of planets was delivered during 
this final phase of planet growth. Similar results have also been reported by O'Brien et al 
(2006), Walsh et al (2011), and Morbidelli et al (2012). This result is particularly important 
in connection with geochemical models of the early differentiation of Earth. As proposed 
by Wood et al. (2008) and Rubie et al. (2011), the oxidation state of Earth continued 
progressively as the Earth was growing. The most successful models of the Earth's core 
formation require that the final 30-40% of Earth's accreted mass to be much more oxidized 
than the material accreted earlier (Rubie et al., 2011). This suggests that the water was 
added together with moderate level of volatile elements during the final 30-40% of Earth's 
accretion and it does not take into account the amount of water that could have been 
adsorbed on dust grains in the early solar nebula (Rubie et al, 2011). Our results (Table 
7) also have direct implication to the recent finding by Marty (2012) who suggested that 
the accretion of volatile elements occurred during the main Earth-forming process instead 
of being a contribution from a late veneer event. 
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3.2.3. Implications of D/H ratio 

As mentioned in the introduction, the primary reason for comets not being considered 
as the main source of Earth's water is the contrast between the measured value of the D/H 
ratio of water in some of the comets with that of SMOW. When comparing these two values, 
it is assumed that the ratio of deuterium to hydrogen in comets is identical to its primordial 
value and the D/H ratio of SMOW is also representative of the time when water was delivered 
to Earth by comets through a late veneer process. Although these two assumptions are not 
entirely valid (radiations and space weathering could have affected the chemical compositions 
of comets during their approach to Earth, the chemistry of SMOW could have changed due 
to the biological processes resulted from the evolution of life, and the collisions of comets with 
Earth could have also affected the chemistry of the material that was transferred to Earth 
from the impactor), our treatment of collisions between two objects (i.e., perfectly inelastic 
with the complete transfer of all material from the impactor to the impactee with no loss of 
water) allows us to evaluate the amount of cometary water that would be necessary to raise 
the D/H ratio of the final objects in the HZ to that of SMOW. In doing so, we consider the 
value of the D/H ratio of a planet, (D/H) P i anet , to be given by a combination of the D/H ratio 
of asteroidal water, (D/H)Ast, and the D/H ratio of the water adsorbed from the primordial 
nebula, (D/H) Neb . The cometary water contribution is not considered at this stage and is 
only included afterwards in order to estimate the amount that is needed to increase the D /H 
ratio of water in the planet to the current value of SMOW. The following equation presents 
a simple relation between the quantities (D/H)Ast, (D/H)N C b and (D/H)pi anet 

K(D/R) Ast + (1 - K)(D/H) Ncb = (D/H) Planet . (5) 

The value of the D/H ratio of asteroids, (D/H)Ast, varies between 1.2 x 10~ 4 to 3.2 x 1CT 4 
(Lecuyer et al. 1998). Adopting the value used in our simulations, (D/H)Ast = 2.2 x 10~ 4 , 
and using the D/H ratios of SMOW and nebula as given in Table 1, Eq. (5) suggests that 
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in order for the asteroidal water to complement the value of the D/H ratio of the nebula 
and raise the D/H ratio of water in the final planets to the current value of SMOW, 64% of 
the amount of the water in these objects has to be delivered by asteroids. In other words, 
a planet that receives 64% of its water from asteroids (and, therefore, the remaining 36% 
from water adsorbed in the solar nebula) will have a D/H ratio similar to that of Earth 
(SMOW). Such a planet will not need any cometary water to raise the value of its D/H 
ratio. Our simulations show that only three planets in model A (simulations 12A-J, 13A-J, 
19A-JS, see Tables 4, 5 and Al) received more than 64% asteroidal water and their D/H 
ratio became higher than that of SMOW. Other interesting cases are, for instance, planets 
produced in simulations 4A-J, 23A-JS, 27A-JS whose water contents, semimajor axes and 
orbital eccentricities are close to those value of Earth (Tables 4, 5 and Al), and they still 
need approximately 2.84%, 11.03% and 20.45% cometary water in order for their D/H ratios 
to reach to that of SMOW. 

The simulations of model B also present interesting results. From the total of 39 planets 
that formed in the HZ, 11 were from these simulations. Table 5 shows the results. Since in 
the simulations of model B, the disk contained only planetary embryos, the cometary water 
necessary to raise the D/H ratio of these planets to that of SMOW varied between 0% and 
45%. As shown by Table 5, in four simulations, the final planets had water contents with 
D/H values larger than that of SMOW (simulations 4B-J, 9B-J, 15B-JS, 17B-JS, see Tables 
6 and Bl). 

An interesting result, when focusing only on the final planets around 1 AU (model A) 
is that the D/H ratio is higher in simulations with only Jupiter compared to those where 
both Jupiter and Saturn were included. This can be related to the strong effect of Saturn, 
in particular the secular resonances, in dynamically exciting the orbits of planetesimals and 
planetary embryos which affects the efficiency of radial mixing from the outer region of the 
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disk to the inner parts. Results of our simulations also indicate that despite a large amount 
of cometary material that is needed in a few cases to increase the D/H ratio of the final 
planets in the HZ to that of SMOW (Tables 4, 5 and 6), the amount of necessary cometary 
water in some cases agrees very well with previous studies [e.g., < 10% as in Morbidelli et 
al. (2000), up to 12% as in Deloule, Robert & Doukhan (1998), and up to 15% as in Owen 
& Bar-Nun, (2000)] giving relevance to locally adsorbed water. 

Although the calculated amount of the cometary water necessary to increase the D/H 
ratio of the final planets in the simulations of models A and B agree with estimates from 
other works, an important questions would be whether models of the dynamical evolution 
of the solar system would be able to accommodate the delivery of this amount of cometary 
material to Earth. As can be seen from the results shown in Tables 4 and 5, adding 0.06O ffi 
of cometary water as estimated by Gomes et al. (2005) or even ~5x 10 _5 M e of cometary 
material as calculated by Morbidelli et al. (2000) to the final planets of these simulations in 
the HZ would not change the value of their D/H ratios significantly. If we also consider the 
loss of water due to the impact of bodies, which according to Marty & Yokochi (2006) could 
be as high as 20%, the situation becomes even worse. An interesting result, however, is that 
as shown by the simulations of model B using our compound model, it is still possible to form 
planets almost only carrying asteroidal water (~ 80%) and with D/H ratios similar to that 
of SMOW. Simulations show that the majority of such objects will obtain their water from 
planetary embryos that originally resided at distances beyond 2.5 AU. The latter implies 
that the initial distribution of water in the protoplanetary disk plays an important role in 
the final D/H ratios of terrestrial planets in the HZ. 
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4. Concluding Remarks 

Assuming that Earth's water had likely more than one source, we studied the contri- 
bution of the water locally adsorbed from the nebula to the water contents and D/H ratio 
of terrestrial planets around the Sun. We carried out 54 numerical simulations considering 
locally adsorbed and asteroidal water, and analyzed the final distribution of mass and water 
content of the planets formed in the Sun's HZ. Results of our simulations indicated that 
large planets with large amounts of water can form in the region between 0.5 and 1.5 AU. 
Results also suggested that our compound model seems to be more efficient in the amount 
and time of the delivery of water over a model in which water is transferred through mere 
collisions of protoplanetary bodies. 

We also studied the D/H ratios of the final planets formed in the Sun's HZ and deter- 
mined the amount of the cometary material that would be required to raise the values of 
their D/H ratios to that of SMOW. Results indicated that using our compound model of 
water distribution and assuming the D/H ratios of comets are representative of their pri- 
mordial values, on average 20% (individually, in some cases, less than 10%) of cometary 
water would be necessary to raise the values of the D/H ratio in the final planets to that 
of SMOW. However, models of the dynamical evolution of the solar system do not seem to 
be able to deliver this amount of cometary material to the HZ of the Sun. For instance, as 
shown by Morbidelli et al. (2000), comets can only provide up to 10% of the Earth's water. 
Our simulations indicate that such a lack of cometary material presents no barrier to the 
formation of terrestrial planets with orbital elements, water contents, and D/H ratios similar 
to those of Earth. Several of our simulations were able to produce Earth-like planets in the 
HZ. 

In carrying out our simulations, we made certain simplifying assumptions. In order to 
avoid complications with breakage and fragmentation of colliding bodies, we assumed that 
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all collisions were perfectly inelastic. We also assumed that in a collision, the mass and the 
amount of the water of the final object would be equal to the sum of the masses and water 
quantities of the impacting bodies. As pointed out by Haghighipour & Raymond (2007), 
this is an assumption that sets an upper limit for the water budget of terrestrial planets 
and ignores the loss of water due to the impacts (Marty & Yokochi 2006) or hydrodynamic 
escape (Matsui & Abe 1986). In a more realistic model, the loss of volatile materials in large 
impacts has to be taken into consideration (Genda & Abe 2005, Canup & Pierazzo 2006). 

Our assumption on the temperature profile of the disk is also one of the limitations of our 
study. The amount of the adsorbed water onto the solar nebula grains and the distribution of 
asteroidal water in the protoplanetary disk strongly depends on the disk temperature profile. 
In our simulations, we followed the model presented by Clark (1998). However, there is a 
rich literature on disk temperature profile and models of water adsorption in accretion disks 
(e.g., Boss 1996, Clark 1998, Sasselov & Lecar, 2000, Fegley 2000, Muralidharan et al. 2008, 
Albarede, 2009). The choice of different temperature profile will clearly result in planetary 
bodies with different amounts of water. The analysis of the connection between the disk 
temperature and the water contents of the final terrestrial planets is outside the scope of 
this paper, and will be presented in future articles. 

Another limitation of our study is its low resolution. Since the speed of an A-body 
integration is proportional to N 2 , in order to keep the time of our simulations in a manageable 
level, we limited our integrations to only a few hundred objects. As such, the results may 
not be able to reveal detailed characteristics of the terrestrial planets of our solar system. 
For instance, our simulations are not able to reproduce the small eccentricities of these 
objects. As shown by Agnor, Canup & Levison (1999) and Chambers (2001), low-resolution 
integrations can produce the main general properties of the final assembly of the planetary 
bodies, however, high resolution simulations such as those by O'Brien et al. (2006) and 
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Raymond et al. (2006, 2007, 2009) are necessary in order to reproduce detailed dynamical 
properties of planets, such as the small eccentricities of Earth and Venus. 

In general, our simulations produced between two and four planets (a < 2AU) on well- 
separated and stable orbits. Venus-size planets are formed in most of our simulations orbiting 
at around 0.5 AU. Several of our simulations have produced a planet inside the habitable 
zone of the system in a timescale of 50-150 Myr, consistent with the expected time of the 
formation of Earth (Jacobsen 2005; Touboul et al 2007). Mercury-size planets did not form 
in our simulations due to our choices for the initial individual masses of planetary embryos 
distributed in the inner part of disk and the position of the inner edge of protoplanetary disk 
(Hansen, 2009). Such difficulties (e.g., the formation of Mars- and Mercury-size objects and 
the architecture of the asteroid belt) have been known to exist in the simulations of terrestrial 
planet formation using similar initial conditions (e.g., Wetherill, 1991; Morbidelli et al., 2000; 
Chambers, 2001; Raymond et al., 2004; 2006; 2007; 2009; O'Brien et al., 2006). Only recently 
efforts have been made to develop models that address these difficulties (Hansen 2009; Walsh 
et al 2011, Izidoro et al. 2013). High- resolution simulations would be necessary to test our 
compound model against these complexities as well. 

As discussed in this paper, no sole source of water provides a satisfactory explanation 
for the origin of Earth's water. The main argument against the water adsorption process 
(and our compound model) as a mechanism for contributing to the Earth's water lies in the 
association of different classes of meteorites to different taxonomic type of asteroids, and 
the connection between their corresponding water contents and their heliocentric distance. 
Carbonaceous chondrites are associated with C-type asteroids, which are primarily beyond 
2.8 AU. Ordinary chondrites with 0.1% of their mass as water are considered fragments of S- 
type asteroids, which are between 2 and 2.5 AU. Enstatite chondrites which are very dry with 
only 0.01% water (Abe et al, 2000) are linked to E-type asteroids at 1.8 AU. All this indicates 
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that the building blocks of Earth around 1 AU must have been extremely dry. However, as 
proposed by Drake (2005), it is uncertain whether the parent bodies of ordinary meteorites 
were indeed anhydrous, or if these meteorites were derived from the metamorphosed outer 
parts of hydrous asteroids. A deep understanding of the relationship between meteorites, 
their parent bodies, and taxonomic type of asteroids is extremely important to improve 
models of origin of Earth's water. We believe that missions such as Dawn and OSIRIS-REx 
will be able to clarify this point. 

Finally, our assumptions on the D /H ratio of comets and that of SMOW present another 
limitation of our study. We assumed that the current D/H ratio of Earth's ocean water is 
representative of its primordial value. However, this quantity could have changed over the 
time due to biological (life and its evolution) and chemical processes. Campins, Swindle & 
Kring (2004) pointed out that the processes involved in planetary accretion such as degassing, 
and the evolution of hydrosphere and atmosphere are complex and may have fractionated 
the chemical and isotopic signature of the source (s) of water. There are also debates on 
whether the estimated bulk D/H ratio and bulk abundance of water on Earth is indeed 
truly known (Abe et al. 2000; Drake and Righter 2002; Smith et al. 2006; Genda & Ikoma 
2008). We assumed that the D/H ratio of SMOW was primordial because of the limited 
knowledge of the evolution of D/H ratio in Earth's water, its relation to its primordial value 
(Williams & Hemley, 2001), and the lack of a sizable statistical sample of the D/H ratio of 
water in other bodies of our solar system. To improve the analysis presented here requires 
wide understanding of the D/H ratio in comets, asteroids and meteorites, a detailed model of 
the variation of D/H ratio during collisions between these bodies, as well as the modeling of 
atmospheric escape and fractionation in Earth's core (Villanueva et al. 2009). Measurements 
of the D/H ratio in a larger sample of comets with HIFI (Heterodyne Instrument for the 
Far Infrared) on Herschel space telescope can also help constraining Solar nebula models 
(Hartogh et al. (2009). 
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Fig. 1. — : Final orbital configurations of 18 simulations of model A considering only Jupiter. The values of the 
disk surface density at 1 AU are shown on the vertical axis. The solar system is also shown for a comparison. 
The size of each body corresponds to its relative physical size, however, it is not to scale on the x-axis. The color 
of each planet represents its water-mass fraction. The eccentricity of each body is represented by a horizontal 
line depicting its variation in its heliocentric distance. 
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Fig. 1. — : Final orbital configurations of 18 simulations of model A considering Jupiter and Saturn. The 
values of the disk surface density at 1 AU are shown on the vertical axis. The solar system is also shown for a 
comparison. The size of each body corresponds to its relative physical size, however, it is not to scale on the 
x-axis. The color of each planet represents its water-mass fraction. The eccentricity of each body is represented 
by a horizontal line depicting its variation in its heliocentric distance. 
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Fig. 2. — : Final configuration of 18 simulations of model B considering only Jupiter (SIM1B-J to Sim9B-J) and 
Jupiter-Saturn (SimlOB-JS to Siml8B-JS), for different values of the disk's surface density at 1 AU. The solar 
system is shown for a comparison. The size of each body corresponds to its relative physical size, however, it is 
not to scale on the x-axis. The color of each planet represents its water-mass fraction. The eccentricity of each 
body is represented by a horizontal line depicting its variation in its heliocentric distance. 
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Fig. 3. — : Initial and final mass of embryos and planetesimals in all the simulations. Upper left: Model A 
with only Jupiter. Upper right: Model A with Jupiter and Saturn. Lower left: Model B with only Jupiter. 



Lower right: Model B with Jupiter and Saturn 
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Fig. 4. — : Final semimajor axis versus the final mass (-M©) of the planets in the simulations of model A. Note 
that each box shows the final results of all 6 simulations corresponding to its surface density (at 1 AU) and 
giant planets configurations as explained in the text. The size of each body corresponds to its relative physical 
size, however, it is not to scale on the s-axis. The color of each planet represents its water content. 
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Fig. 5. — : Final semimajor axis versus the final mass (M®) of the planets in the simulations of model B. Note 
that each box shows the final results of all 3 simulations corresponding to its surface density (at 1 AU) and 
giant planets configurations as explained in the text. The size of each body corresponds to its relative physical 
size, however, it is not to scale on the s-axis. The color of each planet represents its water content. 
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Table 1:: Values of the D/H ratio for comets, nebular gas, and SMOW. 



Body 


D/H Ratio 


Reference 


Halley 


3.16±0.34 x 10" 4 


Eberhardt et al. (1995) 


Hyakutake 


2.9±1.0 x 10~ 4 


Bockelee-Morvan et al. (1998) 


Hale-Bopp 


3.3±0.8 x 10~ 4 


Meier et al. (1998) 


C/2001 Q4 


4.6±1.4 x 10~ 4 


Weaver et al. (2008) 


C/2002 T7 


2.5±0.7 x 10~ 4 


Hutsemekers (2008) 


8P/Tuttle 


4.09±1.45 x 10" 4 


Villanueva et al. (2009) 


153P/Ikeya-Zhang 


< 2.8 ±0.3 x 10~ 4 


Biver et al. (2006) 


C/2004 Q2 


< 2.2 x 10~ 4 


Biver et al. (2005) 


Hartley 2 


1.61 ±0.24 x 10- 4 


Hartog et al. (2011) 


Nebula 


2.1±0.4 x 10~ 5 


Lellouch et al. (2001) 


SMOW 


1.49±0.03 x 10~ 4 


Lecuyer, Gillet & Robert (1998) 



Table 2:: Initial total mass and the amount of mass lost due to different mechanisms in simulations of model 
A (J indicates only Jupiter and JS represents Jupiter-Saturn configuration). From left to right, the columns 
show the simulation's number, surface density at 1 AU, initial total mass of the disk, total ejected mass from 
the disk, the amount of the mass collided with Sun, the amount of the mass collided with Jupiter, and the 
amount of the mass collided with Saturn. The unit of mass is M®. 



Sim 


Si 


Total 


Ejected 


Collisions with 


Collisions with 


Collisions with 






Mass 


Mass 


Sun 


Jupiter 


Saturn 


1A-J 


6 


3.8480 


1.7684 


0.0453 


0.1778 


- 


2A-J 


8 


4.4677 


2.2501 


0.0200 


0.2158 


- 


3A-J 


10 


5.4021 


2.2587 


0.0100 


0.2922 


- 


4A-J 


6 


3.9221 


2.0949 


0.0309 


0.1300 


- 


5A-J 


8 


4.6211 


2.1437 


0.0759 


0.2000 


- 


6A-J 


10 


5.3003 


2.1920 


0.0681 


0.2200 


- 


7A-J 


6 


3.8642 


1.8874 


0.0300 


0.2804 


- 


8A-J 


8 


4.3866 


1.9663 


0.0000 


0.2464 


- 


9A-J 


10 


5.1552 


2.7811 


0.0200 


0.1800 


- 


10A-J 


6 


3.9105 


1.9366 


0.0100 


0.2333 


- 


11A-J 


8 


4.6122 


2.1914 


0.0100 


0.1970 


- 


12A-J 


10 


5.1824 


2.3334 


0.0100 


0.2600 


- 


13A-J 


6 


3.9518 


2.0945 


0.0100 


0.1300 


- 


14A-J 


8 


4.6020 


2.0644 


0.0200 


0.2382 


- 


15A-J 


10 


5.1811 


2.4240 


0.0200 


0.1500 


- 


16A-J 


6 


3.9563 


1.7749 


0.0321 


0.1700 


- 


17A-J 


8 


4.5001 


2.3623 


0.0000 


0.1700 


- 


18A-J 


10 


5.2367 


2.2752 


0.0100 


0.3377 


- 


19A-JS 


6 


3.9516 


2.4910 


0.0300 


0.0200 


0.00 


20A-JS 


8 


4.5495 


2.5844 


0.1576 


0.0652 


0.01 


21A-JS 


10 


5.2066 


2.8992 


0.0980 


0.0900 


0.00 


22A-JS 


6 


3.9444 


2.1014 


0.0426 


0.0500 


0.01 


23A-JS 


8 


4.5990 


2.5392 


0.0100 


0.0300 


0.00 


24A-JS 


10 


5.2309 


2.9351 


0.0100 


0.0600 


0.00 


25A-JS 


6 


4.0377 


2.3283 


0.0982 


0.0600 


0.01 


26A-JS 


8 


4.4341 


2.4578 


0.0400 


0.0700 


0.00 


27A-JS 


10 


5.1971 


2.7734 


0.1220 


0.0900 


0.00 


28A-JS 


6 


3.9163 


1.9476 


0.0100 


0.0721 


0.00 


29A-JS 


8 


4.5052 


2.5096 


0.0200 


0.1654 


0.00 


30A-JS 


10 


5.3477 


3.0595 


0.0200 


0.0700 


0.02 


31A-JS 


6 


3.8367 


2.2518 


0.1190 


0.0400 


0.00 


32A-JS 


8 


4.4984 


2.6605 


0.1018 


0.0500 


0.00 


33A-JS 


10 


5.1142 


2.6800 


0.0757 


0.1000 


0.00 


34A-JS 


6 


3.9050 


2.4512 


0.0839 


0.0516 


0.01 


35A-JS 


8 


4.6589 


2.3875 


0.0200 


0.0600 


0.00 


36A-JS 


10 


5.1265 


2.4865 


0.2034 


0.0700 


0.00 
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Table 3:: Initial total mass and the amount of mass lost due to different mechanisms in simulations of model 
B (J indicates only Jupiter and JS represents Jupiter- Saturn configuration). From left to right, the columns 
show the simulation's number, surface density at 1 AU, initial total mass of the disk, total ejected mass from 
the disk, the amount of the mass collided with Sun, the amount of the mass collided with Jupiter, and the 
amount of the mass collided with Saturn. The unit of mass is M®. 



Sim 


Si 


Total 
Mass 


Ejected 
Mass 


Collisions with 
Sun 


Collisions with 
Jupiter 


Collisions with 
Saturn 


1B-J 


6 


2.8853 


1.3513 


0.0000 


0.0918 




2B-J 


8 


3.9657 


1.8604 


0.0000 


0.0627 




3B-J 


10 


5.0178 


1.9455 


0.1388 


0.2843 




4B-J 


6 


3.0589 


1.4883 


0.0723 


0.0356 




5B-J 


8 


4.0483 


2.1117 


0.0000 


0.0000 




6B-J 


10 


4.8143 


2.8193 


0.0000 


0.0650 




7B-J 


6 


3.1053 


1.3486 


0.0740 


0.1417 




8B-J 


8 


4.0984 


1.5419 


0.0282 


0.1470 




9B-J 


10 


5.0616 


2.4504 


0.0000 


0.1088 




10B-JS 


6 


2.9281 


1.6547 


0.0000 


0.0487 


0.00 


11B-JS 


8 


3.9137 


1.9182 


0.1108 


0.0000 


0.00 


12B-JS 


10 


4.8776 


2.4426 


0.0000 


0.0000 


0.00 


13B-JS 


6 


3.0913 


1.7351 


0.0000 


0.0000 


0.00 


14B-JS 


8 


3.8632 


1.9304 


0.0856 


0.0801 


0.00 


15B-JS 


10 


4.8309 


2.0613 


0.0876 


0.1550 


0.00 


16B-JS 


6 


3.0828 


1.5144 


0.0998 


0.0000 


0.00 


17B-JS 


8 


3.9678 


2.1129 


0.1262 


0.0910 


0.00 


18B-JS 


10 


5.0213 


2.8594 


0.0000 


0.0000 


0.11 
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Table 4:: Final planets inside of HZ in simulations of model A considering only Jupiter. From left to right, 
the columns show simulation's number, semimajor axis, eccentricity, mass (Af ffi ), amount of water (O©), 
percentage of asteroidal mass (> 2.5AU), percentage of asteroidal water (> 2.5AU), percentage of cometary 
water needed to raise the D/H ratio to the current value of SMOW, time (Myr) of delivery of 10©, 20©, 
50©, 100© and 150©. For a comparison the values between Parentheses were obtained using the model of 
water distribution as in Raymond et al. (2004; 2006; 2009). When the values obtained using our model are 
equal to those obtained using Raymond's model, only one entry is shown. 



Sim 


Of 


ef 


Mass 


Water 


%M aat 


%W ast 


%W C 


*del 


£doi 


£doi 


*del 


<del 








(M e ) 


(O®) 










to© 


2O 


50 e 


1OO 


150 e 


1A-J 


1.11 


0.11 


0.76 


13.5 (9.0) 


5.3 


63.1 


(95.1) 


1.6 (-) 


7.5 (20.8) 


12.2 (20.8) 


20.8 (221.4) 


221.4 (x) 


X 


2A-J 


0.96 


0.08 


1.18 


16.4 (9.4) 


3.4 


52.1 


(90.5) 


13.9 (-) 


12.0 (20.2) 


20.2 


20.8 (21.4) 


37.6 (x) 


141.7 (x) 


3A-J 


1.05 


0.09 


1.00 


13.7 (7.0) 


3.0 


46.8 


(91.2) 


18.8 (-) 


8.2 (37.9) 


8.7 (37.9) 


37.9 (110.7) 


68.4 (x) 


X 


4A-J 


1.24 


0.01 


1.08 


20.6 (13.2) 


5.6 


62.1 


(96.6) 


2.8 (-) 


13.3 (29.6) 


20.3 (29.6) 


35.0 (232.2) 


232.2 


232.2 (x) 


5A-J 


1.27 


0.09 


0.79 


15.2 (7.5) 


3.8 


42.1 


(85.1) 


22.7 (-) 


0.2 (10.2) 


0.5 (10.2) 


3.3 (68.9) 


68.9 (x) 


141.8 (x) 


9A-J 


1.02 


0.04 


1.38 


13.5 (5.5) 


1.5 


31.6 


(77.6) 


30.2 (-) 


4.6 (13.2) 


7.4 (13.2) 


13.2 (23.0) 


19.7 (x) 


X 


12A-J 


1.29 


0.11 


0.80 


14.6 (11.3) 


6.3 


72.9 


(94.6) 




22.5 


22.5 


64.8 (68.1) 


68.1 (111.9) 


X 


13A-J 


1.08 


0.02 


0.98 


18.7 (13.2) 


6.2 


68.5 


(96.8) 




2.7 


2.7 


37.4 (46.6) 


46.6 (69.5) 


69.5 (x) 


15A-J 


1.04 


0.09 


0.56 


4.5 (2.4) 


1.8 


47.0 


(87.4) 


18.6 (-) 


9.1 (33.8) 


33.8 


X 


X 


X 


15A-J 


1.29 


0.08 


0.67 


13.0 (7.3) 


4.5 


49.2 


(87.3) 


16.6 (-) 


7.8 


7.8 


8.3 (67.3) 


67.3 (x) 


X 


16A-J 


0.96 


0.04 


0.60 


6.1 (2.7) 


1.7 


34.7 


(78.8) 


28.0 (-) 


56.3 (62.1) 


62.1 


62.1 (x) 


X 


X 


17A-J 


1.12 


0.07 


1.03 


14.1 (7.1) 


2.9 


45.4 


(90.0) 


19.9 (-) 


0.2 (51.8) 


12.8 (51.8) 


51.8 (279.5) 


242.2 (x) 


X 


18A-J 


1.37 


0.17 


1.20 


12.2 (7.1) 


2.5 


52.6 


(90.7) 


13.4 (-) 


33.9 


33.9 


82.2 (146.2) 


146.2 (x) 


X 
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Table 5:: Final planets inside of HZ in simulations of model A considering Jupiter and Saturn. From left 
to right, the columns show simulation's number, semimajor axis, eccentricity, mass (Af©), amount of water 
(O©), percentage of asteroidal mass (> 2.5AU), percentage of asteroidal water (> 2.5AU), percentage of 
cometary water needed to raise the D/H ratio to the current value of SMOW, time (Myr) of delivery of 
10©, 2O0, 50©, 10O© and 15(9©. For a comparison the values between Parentheses were obtained using 
the model of water distribution as in Raymond et al. (2004; 2006; 2009). When the values obtained using 
our model are equal to those obtained using Raymond's model, only one entry is shown. 
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10O® 
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75.5 (227.3) 


X 


X 


21A-JS 
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0.14 


0.85 


9.5 (4.6) 


2.4 


44.9 (93.5) 


20.4 (-) 


8.7 (27.6) 


27.6 


52.5 (x) 


X 


X 


22A-JS 


1.31 


0.02 


0.55 


5.4 (2.3) 


1.8 


39.6 (93.7) 


24.6 (-) 


20.1 (67.6) 


52.7 (67.6) 


101.5 (x) 


X 


X 


23A-JS 


0.95 


0.02 


0.74 


7.8 (4.3) 


2.7 


54.9 (100.0) 


11.0 (-) 


18.3 (25.1) 


25.1 


30.2 (x) 


X 


X 


24A-JS 


0.99 


0.04 


1.47 


12.3 (3.1) 


0.7 


17.3 (68.8) 


38.3 (-) 


6.9 (30.5) 


9.8 (30.5) 


16.9 (x) 


43.2 (x) 


X 


25A-JS 


1.09 


0.03 


0.60 


7.0 (2.4) 


1.7 


30.3 (87.6) 


30.9 (-) 


12.9 (107.5) 


48.1 (107.5) 


107.5 (x) 


x 


X 


26A-JS 


1.27 


0.19 


0.61 


8.1 (4.3) 


3.3 


53.0 (100.0) 


13.0 (-) 


16.1 


16.1 


40.6 (x) 


X 


X 


27A-JS 


1.02 


0.03 


1.07 


14.3 (7.0) 
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44.8 (91.0) 


20.4 (-) 


0.1 (49.3) 


0.2 (49.3) 


40.8 (135.3) 


135.3 (x) 


X 


28A-JS 


0.94 


0.05 


0.44 
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8.5 (-) 


37.2 (63.2) 
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x 


X 


X 
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0.11 


0.77 


6.4 (2.1) 
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28.9 (-) 


13.9 (45.2) 


13.9 (45.2) 


45.2 (x) 


X 


X 


31A-JS 


0.97 


0.12 


0.65 


3.2 (0.1) 


0.0 


0.0 


45.9 


70.9 (x) 


146.4 (x) 


X 


X 


X 


34A-JS 


1.08 


0.02 


0.84 


7.8 (2.3) 


1.2 


27.4 (94.0) 


32.8 (-) 


26.4 (98.2) 


31.6 (98.2) 


98.2 (x) 


X 


X 


35A-JS 


1.22 


0.05 


1.08 


8.9 (0.5) 


0.0 


0.0 


45.9 


3.8 (x) 


35.4 (x) 


52.4 (x) 


X 


X 


36A-JS 


1.01 


0.02 


1.28 


9.0 (2.1) 


0.8 


23.7 (100.0) 


34.9 (-) 


14.3 (17.3) 


14.6 (17.3) 


17.3 (x) 


X 


X 
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Table 6:: Final planets inside of HZ in simulations of model B (J indicates only Jupiter and JS represents 
Jupiter-Saturn configuration). From left to right, the columns show simulation's number, scmimajor axis, 
eccentricity, mass (M®), amount of water (O®), percentage of asteroidal mass (> 2.5AU), percentage of 
asteroidal water (> 2.5AU), percentage of cometary water needed to raise the D/H ratio to the current 
value of SMOW, time (Myr) of delivery of 10®, 20®, 50®, 100® and 150®. For a comparison the values 
between Parentheses were obtained using the model of water distribution as in Raymond et al. (2004; 2006; 
2009). When the values obtained using our model are equal to those obtained using Raymond's model, only 
one entry is shown. 
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0.0 (0.0) 
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X 
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4.3 
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11.7 (-) 
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X 
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41.9 (x) 
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17B-JS 
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Table 7:: Analysis of the amount and time of water-delivery in planets inside of HZ (see Tables 4, 5 and 6). 
From left to right the columns show the disk model (J indicates only Jupiter and JS represents Jupiter-Saturn 
configuration), water distribution model, amount of water delivered, number of planets that received such 
amount of water, mean water-delivery time, range of the water-delivery time, mean of the 60% accretion 
time, range of the 60% accretion time, mean-time of the last giant collision, and the range of the time of the 
last giant collision. 
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Mean-time of (Myr) 


Mean-time (Myr) Mean-time (Myr) of Last 




Distribution Model 
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delivery (range) 


of 60% accretion (range) Giant Impact (range) 


A-J 






All (13) 




61.2 (9.11 - 232.2) 159.6 (48.87 - 291) 




Our Model 














> 20 ffi 


13 


19 (0.5 - 62.1) 








> 50 ffi 


12 


36 (3.3 - 82.2) 
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11 


110 (19.7 - 242.2 ) 
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4 


146 (69.5 - 232.2 ) 






Raymond's Model 














> 20 ffi 


13 


26 (2.7 - 62.1) 








> 50 ffi 


11 


117 (21.4 - 279.5) 
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2 
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All (14) 




44.4 (16.9 - 98.18) 145.6 (30.2 - 268.8) 




Our Model 
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12 
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2 
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_ 
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12 
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> 50 ffi 


2 


181 (135.3 - 227.3) 








> 10O ffi 
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All (5) 
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4 


29 (0.001 - 81.5 ) 
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3 


41 (1.5 -81.5) 








> 10O ffi 


1 


1.5 (1.5 ) 
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1 


1.5 (1.5 ) 






Raymond's Model 














> 20 ffi 


2 


42 (1.5 - 81.5) 
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2 


42 (1.5 - 81.5) 








> lOOe 


1 


1.5 (1.5 ) 
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1.5 (1.5 ) 




B-JS 






All (6) 




47.17 (26.48 - 71.57) 105.6 (47.9 - 159.6) 




Our Model 














> 20 ffi 


5 


36 (0.6 - 71.6 ) 








> 5O 


3 


23.3 (12.3 -31.2) 
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3 


23.3 (12.3 -31.2) 
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3 


63 (12.3 - 146.6 ) 






Raymond's Model 














> 20 ffi 


3 


23.3 (12.3 -31.2) 








> 50 ffi 


3 


23.3 (12.3 -31.2) 
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3 


23.3 (12.3 -31.2) 








> I50 e 


2 


22 (21.75 ) 





Final configuration of the simulations of model A 
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Table Al:: Final results of the simulations of model A (J indicates only Jupiter and JS represents Jupiter- 
Saturn configuration). From left to right, the columns show simulation's number, surface density at 1 AU, 
semimajor axis, eccentricity mass (M®), amount of water (O®), percentage of asteroidal mass (> 2.5 AU), 
percentage of asteroidal water (> 2.5 AU), final D/H ratio, and the time (Myr) of the delivery of lO®, 20®, 
50 e , 10O© and 150®. 
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0.02 
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2.0 
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1.4ole-U4 


4.8 
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X 


X 


X 


1 A T 

1A-J 


6 
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0.11 


0.76 


13.5 


5.3 


63.1 
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12.2 


20.8 


221.4 


X 
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D 
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2A-J 


8 
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0.16 
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0.0 


0.0 
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X 


X 


X 


X 


X 
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8 


0.96 


0.08 
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12.0 
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20.8 


37.6 
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ZA-J 
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2.19 


0.13 


0.37 
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0.2 


0.2 


2.5 
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X 
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10 


0.63 


0.04 


0.75 
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1.972e-04 


73.4 
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X 


X 
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3.0 
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X 
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X 


X 


X 
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6 
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0.01 
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1.446e-04 


13.3 


20.3 


35.0 


232.2 


232.2 


4 A- J 


6 


2.50 


0.18 


0.09 


5.6 


21.3 


76.3 


1.728e-04 


1.4 


2.2 


2.6 


X 


X 


4 A- J 


6 


2.81 


0.16 
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0.2 


0.0 


0.0 
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X 


X 


X 


X 


X 


5A-J 


8 
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0.5 


3.3 


68.9 


141.8 


5A-J 


8 


2.06 
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0.0 


0.0 
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X 


X 
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10 


0.72 


0.03 


1.67 


24.4 


3.6 


52.5 


1.256e-04 


4.8 


9.9 


9.9 


9.9 


16.9 


6 A- J 


10 


1.65 


0.11 


0.93 


10.6 


2.2 


40.4 


1.015e-04 


6.3 


6.3 


12.3 


54.8 


X 


6 A- J 


10 


2.76 


0.16 


0.13 


5.1 


15.5 


82.9 


1 OTA f\ A 

1.859c- 04 


5.4 


5.4 


54.1 


X 


X 
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6A-J 


10 


2.92 


0.15 


0.04 


0.2 


0.0 


0.0 


2.100e-05 


X 


X 


X 


X 


X 


7 A- J 


6 


0.65 


0.01 


0.63 


10.2 


6.4 


83.5 


1.872e-04 


16.0 


16.0 


125.7 


272.4 


X 


7 A- J 


6 


1.38 


0.03 


0.66 


16.0 


7.6 


66.7 


1.537e-04 


4.6 


4.6 


31.9 


101.5 


230.5 


7 A- J 


6 


2.25 


0.19 


0.06 


4.5 


31.0 


93.8 


2. 077c- 04 


3.5 


3.5 


X 


X 


X 


7 A- J 


6 


2.47 


0.34 


0.06 


0.3 


0.0 


0.0 


2.100e-05 


X 


X 


X 


X 


X 


7A-J 


6 


2.51 


0.08 


0.06 


0.7 


0.0 


0.0 


2.100e-05 


X 


X 


X 


X 


X 


7A-J 


6 


2.79 


0.12 


0.07 


7.1 


44.4 


90.1 


2.003c-04 


2.5 


2.5 


17.0 


X 


X 


7A-J 


6 


2.94 


0.17 


0.04 


2.5 


25.6 


86.8 


1.937e-04 


10.9 


10.9 


X 


X 


X 


8A-J 


8 


0.85 


0.05 


1.52 


24.4 


5.3 


69.8 


1.600c-04 


0.0 


12.1 


38.5 


112.8 


150.5 


8A-J 


8 


1.67 


0.08 


0.27 


9.7 


11.1 


66.2 


1.527c-04 


0.4 


0.4 


149.5 


X 


X 


8A-J 


8 


2.82 


0.13 


0.14 


7.1 


22.0 


89.7 


1.996e-04 


140.9 


140.9 


197.0 


X 


X 


8A-J 


8 


2.97 


0.32 


0.04 


0.5 


0.0 


0.0 


2.100e-05 


X 


X 


X 


X 


X 


9A-J 


10 


0.50 


0.25 


0.42 


0.2 


0.0 


0.0 


2.100e-05 


X 


X 


X 


X 


X 


9A-J 


10 


1.02 


0.04 


1.38 


13.5 


1.5 


31.6 


8.382e-05 


4.6 


7.4 


13.2 


19.7 


X 


9A-J 


10 


1.94 


0.21 


0.16 


1.3 


0.0 


0.0 


2.100e-05 


37.9 


X 


X 


X 


X 


9A-J 


10 


2.73 


0.05 


0.23 


4.0 


4.3 


53.3 


1.270c-04 


0.1 


4.4 


X 


X 


X 


10A-J 


6 


0.57 


0.28 


0.17 


0.1 


0.0 


0.0 


2.100e-05 


X 


X 


X 


X 


X 


10A-J 


6 


0.87 


0.02 


0.74 


16.8 


8.2 


76.2 


1.725c-04 


11.8 


11.8 


33.0 


77.1 


214.1 


10A-J 


6 


1.65 


0.02 


0.61 


7.2 


1.6 


29.6 
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14.4 
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71.4 


X 


X 


11A-J 


8 
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0.9 
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58.8 


X 


X 
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8 
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0.0 


0.0 
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X 


X 


X 


X 


X 
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8 
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0.11 
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11.7 
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X 


X 
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8 


1.98 


0.11 


0.46 
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2.2 


28.4 
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0.3 


7.4 


7.4 


X 


X 


11A-J 


8 


2.79 


0.09 


0.05 


0.5 


0.0 


0.0 


2.100e-05 


X 


X 


X 


X 
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Table Al:: Final results of the simulations of model A (J indicates only Jupiter and JS represents Jupiter- 
Saturn configuration). From left to right, the columns show simulation's number, surface density at 1 AU, 
semimajor axis, eccentricity, mass (M ffi ), amount of water (O©), percentage of asteroidal mass (> 2.5AU), 
percentage of asteroidal water (> 2.5 AU), final D/H ratio, and the time (Myr) of the delivery of 10®, 20 ffi , 
50 e , 10O e and 150®. 
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Table Al:: Final results of the simulations of model A (J indicates only Jupiter and JS represents Jupiter- 
Saturn configuration). From left to right, the columns show simulation's number, surface density at 1 AU, 
semimajor axis, eccentricity, mass (M ffi ), amount of water (O©), percentage of asteroidal mass (> 2.5AU), 
percentage of asteroidal water (> 2.5 AU), final D/H ratio, and the time (Myr) of the delivery of 10®, 20 ffi , 
50 e , 10O e and 150®. 
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Table Al:: Final results of the simulations of model A (J indicates only Jupiter and JS represents Jupiter- 
Saturn configuration). From left to right, the columns show simulation's number, surface density at 1 AU, 
semimajor axis, eccentricity, mass (M ffl ), amount of water (O©), percentage of asteroidal mass (> 2.5AU), 



percentage of asteroidal water (> 2.5AU), final D/H ratio, and the time (Myr) of the delivery of 10®, 20®, 
50 e , 10O e and 150®. 
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B. Final configuration of the simulations of model B 
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Table Bl:: Final results of the simulations of model B (J indicates only Jupiter and JS represents Jupiter- 
Saturn configuration). From left to right, the columns show simulation's number, surface density at 1 AU, 
semimajor axis, eccentricity, mass (M ffi ), amount of water (O©), percentage of asteroidal mass (> 2.5AU), 
percentage of asteroidal water (> 2.5 AU), final D/H ratio, and the time (Myr) of the delivery of 10®, 20 ffi , 
50 e , 10O e and 150®. 
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9B-J 


10 


2.14 


0.19 


0.65 


4.2 


0.0 


0.0 


2.100c-05 


15.1 


30.9 


X 


X 


X 


10B-JS 


6 


0.33 


0.15 


0.30 


2.0 


0.0 


0.0 


2.100e-05 


6.9 


X 


X 


X 


X 


10B-JS 


6 


0.84 


0.10 


0.72 


2.7 


0.0 


0.0 


2.100e-05 


58.6 


150.9 


X 


X 


X 


10B-JS 


6 


1.71 


0.30 


0.04 


0.2 


0.0 


0.0 


2.100e-05 


X 


X 


X 


X 


X 


10B-JS 


6 


3.15 


0.14 


0.04 


9.1 


100.0 


100.0 


2.200c-04 


X 


X 


X 


X 


X 


11B-JS 


8 


0.58 


0.10 


0.54 


0.8 


0.0 


0.0 


2.100e-05 


X 


X 


X 


X 


X 


11B-JS 


8 


0.93 


0.18 


0.34 


1.3 


0.0 


0.0 


2.100e-05 


41.5 


X 


X 


X 


X 


11B-JS 


8 


1.43 


0.11 


0.78 


5.9 


0.0 


0.0 


2.100c-05 


2.4 


30.1 


85.3 


X 


X 


11B-JS 


8 


2.18 


0.11 


0.12 


25.2 


100.0 


100.0 


2.200c-04 


0.9 


0.9 


0.9 


0.9 


0.9 


12B-JS 


10 


0.94 


0.05 


2.26 


38.1 


4.3 


54.3 


1.290c-04 


7.5 


17.7 


31.2 


31.2 


31.2 


12B-JS 


10 


1.55 


0.13 


0.27 


1.0 


0.0 


0.0 


2.100e-05 


5.6 


X 


X 


X 


X 


13B-JS 


6 


0.58 


0.17 


0.39 


0.7 


0.0 


0.0 


2.100e-05 


X 


X 


X 


X 


X 


13B-JS 


6 


0.98 


0.10 


0.63 


3.6 


0.0 


0.0 


2.100e-05 


30.5 


71.6 


X 


X 


X 


13B-JS 


6 


2.73 


0.02 


0.04 


9.1 


100.0 


100.0 


2.200e-04 


X 


X 


X 


X 


X 



-56 - 



Table Bl:: Final results of the simulations of model B (J indicates only Jupiter and JS represents Jupiter- 
Saturn configuration). From left to right, the columns show simulation's number, surface density at 1 AU, 
semimajor axis, eccentricity, mass (M®), amount of water (O®), percentage of asteroidal mass (> 2.5AU), 



percentage of asteroidal water (> 2.5 AU), final D/H ratio, and the time (Myr) of the delivery of 10®, 20®, 
50®, 100® and 150®. 



Sim 


Si 


at 


ef 


Mass (M ) 


Water (O e ) 


%M ast 




D/H 


tdel 


tdel 


tdel 


*del 


*dcl 






















20 ffi 


50 ffi 


l0O ffi 


150® 


14B-JS 


8 


0.83 


0.09 


1.56 


8.0 


0.0 


0.0 


2.100e-05 


8.1 


26.0 


47.4 


X 


X 


14B-JS 


8 


2.74 


0.15 


0.06 


12.7 


100.0 


100.0 


2.200c-04 


X 


X 


X 


X 


X 


15B-JS 


10 


0.64 


0.10 


1.13 


6.4 


0.0 


0.0 


2.100e-05 


2.2 


19.9 


55.3 


X 


X 


15B-JS 


10 


1.20 


0.27 


1.19 


25.2 


7.0 


69.9 


1.600c-04 


0.0 


0.6 


12.3 


12.3 


12.3 


15B-JS 


10 


3.04 


0.17 


0.09 


19.7 


100.0 


100.0 


2.200c-04 


X 


X 


X 


X 


X 


16B-JS 


6 


0.68 


0.22 


0.64 


18.7 


12.1 


87.6 


1.952c-04 


20.7 


75.4 


75.4 


75.4 


75.4 


16B-JS 


6 


1.26 


0.18 


0.49 


2.3 


0.0 


0.0 


2.100c-05 


41.9 


61.2 


X 


X 


X 


16B-JS 


6 


3.28 


0.17 


0.09 


18.9 


100.0 


100.0 


2.200c-04 


0.0 


0.0 


0.0 


0.0 


0.0 


17B-JS 


8 


0.62 


0.05 


0.69 


3.5 


0.0 


0.0 


2.100e-05 


71.8 


71.8 


X 


X 


X 


17B-JS 


8 


0.97 


0.05 


0.60 


15.4 


9.7 


80.7 


1.815c-04 


26.5 


26.5 


26.5 


26.5 


146.6 


17B-JS 


8 


1.59 


0.23 


0.25 


1.5 


0.0 


0.0 


2.100e-05 


17.1 


X 


X 


X 


X 


18B-JS 


10 


0.76 


0.10 


1.76 


42.4 


9.8 


86.3 


1.928c-04 


12.8 


33.2 


55.8 


55.8 


55.8 


18B-JS 


10 


2.51 


0.08 


0.42 


60.9 


66.1 


96.5 


2.131e-04 


0.7 


0.7 


0.9 


0.9 


0.9 



